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Abstract 

The transverse momentum dependence of Hanbury-Brown/Twiss (HBT) in- 
terferometry radii for 2-body correlation functions provides experimental ac- 
cess to the collective dynamics in heavy- ion collisions. We present an analyti- 
cal approximation scheme for these HBT radii which combines the recently 
derived model-independent expressions with an approximate determination 
of the saddle point of the emission function. The method is illustrated for a 
longitudinally boost-invariant hydrodynamical model of a heavy ion collision 
with freeze-out on a sharp hypersurface. The analytical approximation con- 
verges rapidly to the width of the numerically computed correlation function 
and reproduces correctly its exact transverse momentum dependence. How- 
ever, higher order corrections within our approximation scheme are essential, 
and the previously published lowest order results with simple m± scaling be- 
haviour are quantitatively and qualitatively unreliable. 
PACS numbers: 25.75.-hr, 07.60.1y, 52.60. -hh 
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I. INTRODUCTION 



While the total energy involved in a heavy ion collision can be measured directly by parti- 
cle calorimeters, an analogous direct measurement of the locally reached energy density does 
not exist. An indirect determination of the size of the interaction region is possible through 
Hanbury-Brown/Twiss (HBT) intensity interferometry [|T|. However, the interpretation of 
the measured correlations between the produced particles is in general model-dependent, 
and a considerable amount of theoretical effort has recently been spent on the question to 
what extent this intrinsic model dependence can be reduced by a refined analysis [0-[7|. 

Compared to HBT interferometry on stars, the situation in heavy-ion collisions is com- 
plicated by the finite lifetime and the strong dynamical evolution of the particle emitting 
source |]8|-[TI|. As a result the 2-body correlation function is in general characterized by 
different width parameters ( "HBT radii" ) if the relative momentum between the two identi- 
cal particles points in different directions 0. Furthermore, the dynamical expansion of the 
source leads to correlations between the momenta of the emitted particles and their emission 
point (so-called x — K correlations) which in turn generate a characteristic dependence of 
the HBT radii on the total pair momentum 0. 

In this paper we discuss this last issue in some detail. Our work was motivated by 
the simple scaling laws for the HBT radii as a function of the transverse mass of the 
particle pair which were proposed in Refs. ||10|,|3|,|6|J7|JI3|1 . These predictions are based on 
simple models for the emission function, and it is natural to ask which of the predicted 
features are independent of the model and which are not. For example, the model of Ref. 
makes definite assumptions about the shape of the longitudinal and transverse expansion 
flow profiles whose influence on the result for the correlation function, in particular on the 
functional form of the m^-dependence of the correlation radii, is not known. The work 
reported here provides an answer to this question within a certain class of models. More 
importantly, however, we discovered on our way that the simple scaling laws of Refs. [p!0| ,^|6|,|7| 
are based on approximations which are quantitatively unreliable and in some typical cases 
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give even qualitatively misleading results. We will show, for example, that a fit of the exact 
correlation function with the simple analytical expression from Ref. ||T0|] for the longitudinal 
HBT radius Ri gives an estimate for the decoupling time which is too large by a factor ~ 2. 
(This was recently also pointed out in |13|.) The effect of the finite duration of particle 



emission on the difference between the HBT radii in "out" and "side" directions, — R^ 
which, as we will show, is also present in models with sharp freeze-out along a proper-time 
hypersurface, is completely missed by the approximations used in Ref. 0, giving rise to a 
qualitatively wrong m± dependence of the "out" -correlator. We conclude from our findings 
that a quantitative analysis of the m_L-dependence of HBT radii necessitates a numerical 
evaluation of the theoretical expressions for the correlation radii, or requires at least a much 
more sophisticated analytical approximation scheme than so far employed. 

We begin by shortly reviewing the general approach and defining our notation. We 
restrict our discussion to two-particle correlations. We start from the usual abstraction of 
the collision region as an assembly of classical boson emitting sources in a certain space-time 
region [^IT^. Their Wigner transform [p!^|-[T7| defines the emission function S{x,p) which 
describes the probability that a particle with on-shell momentum p {p^ = m?) is emitted 
from the space-time point x. The emission function determines the single particle momentum 
spectrum -Pi(p) = EpdN/d^p = J d'^x S{x,p) as well as the HBT two-particle correlation 
function C(K, q) where K = | (pi -|- P2) = (K_|_, K^) is the average momentum of the two 
particles and q = Pi — P2 their relative momentum. In the plane wave approximation 
the latter is given in terms of the 1- and 2-particle distributions -Pi(p), P2(Pi5P2) and the 
average number (A^) {{N{N — 1))) of particles (particle pairs) produced in the reaction as 



111151,1171 



rrn ^ W i^2(pi,P2) \Jd'xSix,K)e^'^-f 

^Pi'P^J (iV(iV-l))Pi(pi)Pi(p2) ^ Jd^xSix^p,) Jd^xSix^p^)- ^ ■ ' 

In the numerator of the last expression we have introduced off-shell momentum 4-vectors 
K and q for the total and relative momentum of the particles in the pair by defining = 



h{Ei + E2) and = Ei — E2 where Ei = J w? + 'p} . A popular approximation whose 



accuracy was studied quantitatively in consists of putting K on-shell, ~ Ek 
^/m? + K^, and setting pi=p2=K in the denominator: 



\Jd'^xS{x,K) e*«-^ 
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C(K,q)~l + ^^ '-^-^ (1.2) 

The aim of HBT interferometry is to obtain from the measured correlation function C(K, q) 
information about S{x, K) which in turn should characterize the size of the interaction 
region. One usually compares C(K,q) to experimental data by taking recourse to the 
Gaussian approximation 0] 

C(K, q) ^ 1 + exp [-qlRl{K) - qlRl{K) - qfR^K) - 2qiq,Rl{K)] , (1.3) 

which is valid for azimuthally symmetric sources. Here qi, qo, qs denote the spatial compo- 
nents of q in the beam direction ("longitudinal" or 2;-direction) , parallel to the transverse 
components of K ("out" or x-direction), and in the remaining third cartesian direction 
("side" or y-direction), respectively. In practice the last term in ( |1.3| ) which mixes the com- 
ponents qo and qi has usually been neglected because its existence was only recently pointed 
out and confirmed experimentally |]TB|. 



The main purpose of the present work is to investigate the transverse momentum de- 
pendence of the HBT correlation radii i?^(K) in ( |1.3|) . As already mentioned it reflects the 
X — /^-correlations in the emission function S{x, K) generated by collective expansion of the 
source, information that can not be obtained from K-averaged correlation radii. Our analy- 
sis will be based on a specific model for the emission function, but the method is general and 
can later be combined with a comprehensive investigation of the model dependence of 2- 
particle correlations. In the present paper we will compare two different methods to calculate 
the Rf(K): the first consists of fitting the numerically determined HBT correlation function 
C(K, q) to the form ( |1.3| ) in such a way that its half width is correctly reproduced, while in 
the second approach we will use the following model-independent expressions [0,^,^]: 

R's = {y\ 

Rl = {{x-p^tf)-{{x-(3^t))\ 



Rl = ii^ - P±t)iz - ht)) - {{x - (3^t)){{z - (3Lt)) , (1.4) 

which are guaranteed to yield the correct Gaussian curvature of the correlation function 
C(K,q) near q = In Eqs. (0) we defined A = 2Ki/(Ei + E2) ^ Ki/Ex and 

introduced the notation 



For Gaussian sources, the set of equations ( |1.4|) can be alternatively derived by expanding 



(|1 ■ 1|) for small relative momenta q up to second order and re-exponentiating [§,|T31; or by 
making a Gaussian saddle point approximation of the emission function around its maximum 
0. A more general justification of Eqs. ( [1.4|) , which (as we will show) remain rather accurate 
even for non-Gaussian sources, is based on the following parametrization of the emission 
function: 

where the saddle point x(K) is defined as the point where all first derivatives of S{x, K) 
vanish. If one identifies B^^, with the tensor of second derivatives of In ^(a;, i^'), the first 
term on the r.h.s. of Eq. ( |1.6| ) amounts to the Gaussian saddle point approximation around 
X which was used in [§-0]. For non-Gaussian forms of S{x,K), the typical corrections to 
Eqs. ( |1.4D , ( |1.5| ) from 6S{x,K) can be minimized by instead defining B^,^(K) in terms of 
the variance of S{x, K), 

(B^^)^^ = {x^x^) - {x^){x^) . (1.7) 

which measures the width of S{x,K). With this definition the Fourier transform ( |1.2| ) of 
the first term in ( p..6| ) can be done analytically and leads directly to the expressions (|1.3| ) 



and ( |1.4| ), with the expectation value ( |1.5| ) defined in terms of the full emission function 
S{x,K). We will see that, in contrast to other approximation schemes (see for example 
[^,1^), the determination of the HBT radii via (|1.4| ) yields an accurate representation of the 



width of the exact correlation function, independent of the vahdity of the Gaussian saddle 
point approximation for S{x,K). 

Except for the few special cases for which the four-dimensional Fourier transform of 
S{x, K) is known exactly, an analytical investigation of equation { \L . 1| ) or ( |1.4[ ) must involve 
a suitable approximation scheme. The existing analytical calculations [^,^,|^,|^ of (the K- 
dependence of) the HBT radii for particular emission functions S{x, K) cover certain limiting 
cases of parameter space only, and their validity has not been checked numerically. In the 
present work, we investigate the complete K-dependence of the HBT radii by both analytical 
and numerical methods without any restriction to limiting cases. 

Our paper is organized as follows: In section we review a recently introduced simple 
model for the emission function which we will also study here. It may not be the 
most realistic model for heavy-ion collisions, but it possesses a number of essential physical 
features and allows for a controlled investigation of our analytical approximation scheme and 
a direct comparison with previously published approximations. In section ^TT| we develop 
a new analytical approximation scheme for the HBT radii ( |1.4D of this model. In Section 
[IV| we provide some intermediate results for the numerical evaluation of the correlation 
function. In Section |V| we compare the analytical and numerical values for the HBT radii. 
We find that the analytical approach works very well for a linear transverse flow profile, if 
the approximation scheme (which involves an expansion in powers of 7t^«_L/T, 7^ being the 
Lorentz factor associated with the transverse flow) is carried to third non-leading order. For 
the "out" and "longitudinal" radii the leading terms in this approximation scheme are found 
to be insufficient. For parabolic transverse flow profiles our approximation scheme and, a 
fortiori, all previously suggested simpler approximations are found to fail completely. The 
physical consequences of our findings are discussed in Section |V|. 
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II. A SIMPLE MODEL FOR THE EMISSION FUNCTION 



For ease of comparison with published results in the literature, we consider the emission 
function 

= :^:!lMe-^VTg-.V(2iJ^) , (2.1) 

[2'Ky 

where the Boltzmann factor q-^-'^I'^ reflects the assumed local thermal equilibration of a 
source with local temperature T{x) moving with hydrodynamic 4- velocity u^{x). We will 
take T to be constant. We assume sharp freeze-out of the particles from the thermalized 
fluid along a hypersurface and the 4-vector n^(x) = JY;d^(T^{x') S^'^\x — x') denotes 

the normal-pointing freeze-out hypersurface element. The factor exp [— r^/(2i?^)] imposes a 
Gaussian transverse density profile with geometric radius i?; the exponent can be interpreted 
in terms of an effective chemical potential fi{r)= — Tr^/(2i?^) which depends only on the 
transverse coordinate r=^x^ + y^. 

For the flow velocity profile u^[x) we assume longitudinal boost invariance by setting 
VL=z/t, i.e. identifying the flow rapidity ?7fiow=| ln[(l + vl)/{^ — vl)] with the space-time 
rapidity ?7=| ln[(t + z) /{t — z)]. We can thus parametrize u'^{x) in the form 

X y 

Ui,{x) = {chrjchrit, shr^^— , sh.r]t-, sh i] ch r]t) , (2.2) 

r r 

where rit{r) is the transverse flow rapidity. The momentum K is parametrized in 



the usual way through the transverse mass m±=^m'^ + Kj_ and the rapidity Y as 
K^={m±ch Y, K±, 0, m^sh Y). 

Assuming sharp freeze-out at constant proper time tq, the freeze-out hypersurface is 
parametrized as S(a;) = (rochr7, x, y, TQshi]), resulting in K-n{x)=m±ch {rj — Y) 6{t — tq). The 
emission function can thus be rewritten as 

S{X,K) = "^-LCh(r7-y) ^^[m|ch(»;-y)chr?t-JCi(a:/r)shr?t]/T^-rV(2fl^) (2.3) 

' (27r)^ 

The 5-function renders the r-integration in d'^x=T dr drj dx dy trivial, and only the inte- 
grations over space-time rapidity t] and the transverse coordinates x, y remain to be done. 



Exploiting the boost invariance and infinite longitudinal extension of our source we can sim- 
plify the structure of Eqs. ( p..4[ ) by going to the longitudinally comoving system (LCMS). In 
this frame Y=0=i3l- 



III. ANALYTICAL EVALUATION OF THE HBT RADII 

In this section we will discuss the model-independent expressions ( |1.4|) for the HBT 
radii. Our analysis of these expressions proceeds in two steps: First, the ?7-integration is 
done analytically. We then find the approximate saddle point of the resulting integrand in 
the X- and |/-directions and carry out those integrations by saddle point integration. 



A. Integration over rj 

To begin with we note that by the symmetry of our infinite boost-invariant emission 
function ( |2.3| ) the expectation values of all odd functions of r] vanish. Therefore, in this 
model the cross term Rio is zero 0. 

Writing t=Tch.ri, z=Tsh.ri and including the ch rj prefactor in Eq. ( |2.3|) (recall that 
we use the LCMS frame where i3l=Y=0), the r^-integrations in the expectation val- 
ues (|1.4| ) can be done analytically and expressed in terms of modified Bessel functions 
K^(a) = J^dr]e-^'''^'^ch{ur]). Defining 



d77sh'"r7ch'^+Se"°^''^'^V / dr]chr]e-''^'^''^'^ (3.1) 

-oo / J — CO 



with 



a{r) = —chrjtir) , (3.2) 
the expressions (|1.4|) for the HBT radii can be written as sums of terms of the general form 

where f{x,y) is a simple polynomial in x or ?/ and 
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is the transverse weight function for the (. . .)* average in (|3.3|). The functions C"" needed 



m 



L4D are given exphcitly by 



a Ki{a) 
G''\r) = 1 



02/ N _ 1 , ^2(0) 



a Ki(a) 



aKi[a) 

with a = a{r) from ( |3.2| ). With these definitions the HBT radii ( |1.4| ) take the form 

?2 2 /a^20 



?2 /_2\ /_\2 



-2/5^ro((a;G°^(r))^-(a;),(G°^(r))J 



Rl = . (3.6) 

For the model ( |2.1| ) these expressions are exact. 

In some previous papers P,[7|,pi| the 77-integration was not done exactly, but by sad- 
dle point integration around f/ = 0. To this end one expands the exponent of ( |2.3| ) as 
chi] ~ 1 + rf/^, and replaces all chi] prefactors by their saddle point value chf/ = 1. 
This amounts to replacing the modified Bessel functions by their leading term for large a, 



K„{a) ~ yn/{2a)e~°'. It is easy to see that in this approximation the ro-dependent terms 
in Eq. (|3.6|) for vanish exactly. One would thus conclude that for sharp freeze-out at 
r = To there is no influence of the source lifetime on R^ and on the difference R^ — R^. 

The exact expressions ( p.6| ) show, however, that this is not true: R^ — R^ is sensitive to 
the time structure of the source. The physical mechanism behind this is quite interesting: 
since particles from different points on the freeze-out surface contribute to the correlation 
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function as long as they are separated in z-direction by less than the longitudinal region of 
homogeneity Ri, the correlation function indeed probes a finite range of coordinate times t 



along the proper-time hyperbola t=tq: At ~ y r^+Rf /A — tq (assuming that only half of 
the longitudinal region of homogeneity contributes effectively). This shows up in the HBT 
radius for the out-direction as a contribution from a finite duration of the emission process 
and is formally reflected by the last two lines in the expression (|3.6| ) for R"^. We will see the 
practical consequences of this in Sections and 



B. Transverse saddle point integration 

Since in most realistic situations we can assume a(r) > ^ > 1 (\ow-K± photons from 
a high-temperature source being the exception), it makes sense to use in the transverse 
integrals the asymptotic expansion of the Bessel functions for large arguments: 

Clearly, how many orders p must be taken into account to obtain a good approximation for 
(|3.3|) depends on both the order v and the lower limit for a, as reflected by the upper bound 
20| 



(2a)Pp!|r(z/-p+|)| ' ^-"^ 2 



|i?(z/,p,a)|< ,^_,J„^Y ' , iM ^ P>^-o, (3- 



for the remainder term in Eq. (|3.7| ). 

Substituting the expansion (|3.7| ) back into (|3.6|) , we are left with sums of integrands of 
the type f{x,y) Gn{r) where 

Gnir) = (Chr/t)-"-^ e-("'±chvt-K^{x/r)shvt)/T ^-ry{2R^) _ ^g^g) 

The integer n labels the order in the expansion (|3.7| ). The resulting transverse integrals 
can be done by a saddle point approximation of Gn{r). To obtain accurate results it is, 
however, important to find a good approximation for the exact saddle point x„ of Gn{r) 
which cannot be obtained analytically except for simple limiting cases P]. In previous 
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studies a Gaussian approximation of tlie first exponential factor in ( |3.9| ) around its 

own saddle point was used; the resulting product of two Gaussian factors was integrated out 
analytically, without however including the effects from the (chrit)'^"''^^^ prefactor in (|3.9|). 



Unfortunately, we found that in most cases the resulting analytical expressions for the HBT 
radii do not exhibit the correct m±_ dependence of the exact results obtained by evaluating 
the expressions ( p. 61) numerically. In Appendix ^ we derive an analytical approximation 
scheme for obtaining a better estimate of the true saddle point which yields much more 
accurate results. We derive there the following Gaussian approximation for the transverse 
weight function G'n(r): 



G„(r) ^ Cn exp 



with effective Gaussian size parameters 



2Rl{n) 2Rl{n) 



(3.10) 



1 1 1 

^ + 



i?2(^) A2(a;„) ' 

Here Xn is the x-component of the transverse saddle point r„ which should be determined 
from the condition 

((-^?)-^«-^-'^«)|-l^- (-2) 

Barred quantities with an index n denote values at the saddle point, and primes indicate 
a;-derivatives. Since this equation cannot be solved analytically for x„, one replaces in 
Eqs. ( |3.11|) the true saddle point x„ by an analytically obtained approximate value x„ and 
calculates the 'lengths of homogeneity' A^(x„) and A„(x„) by evaluating the expressions 
(|A7| , [A8|) at that point. The correct saddle point x„ is then related to the approximate value 
Xn by Eqs. (|A12| , |A13D , while the constant C„ is determined by Eq. ( [A15| ). A good choice 



for the approximate saddle point x„ is also derived in Appendix ^ and given analytically 
by Eqs. P16D , (^ and (|ATO| ). 



Please note that according to Eqs. ( |A7| , [A8D the homogeneity lengths in the 'side' and 
'out' directions. As and Ao, are in general different; this was also recently found by Akkelin 
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and Sinyukov Only for a linear transverse rapidity profile with a small slope (the 
nonrelativistic transverse expansion scenario studied in Refs. [§,0), where the transverse 
saddle point can be taken as x„ ~ 0, do the two expressions agree. We will see in Section |V] 
that in some situations this difference can be drastic, causing sizeable differences between 
Ro and Rg which have nothing to do with a finite duration of the emission process. 

Using the expansion (|3.7|) , the HBT radii (|3.6| ) can be rewritten in terms of integrals 
over the auxiliary functions Gn{r), Eqs. (|3.9| / |3.1(3|) . Both the numerator and denominator 



of ( p.3|) reduce to integrals of the type 

p / J" /. 

^Cn^—j J dxdy f{x,y)Gn{r) , (3.13) 



n=0 



with certain coefficients c„ which depend on the particular HBT radius under consideration 
(see Appendix Q). For the polynomials f{x,y) occuring in ( p.6|) the necessary integrals are 



easily evaluated in terms of the size parameters ( |3.11 



J dxdyGn{r) =2nCnRo{n) Rsiji) , 
j dx dy y^ ^^(r) = 2nGn Rain) Rl{n) , 
J dx dy X Gn{r) = 2iTGn Roin) Rs{n)xn, 

j dxdyx^Gnir) =27TGnRo{n)Rs{n) [Rl{n) + xl) . (3.14) 

Before combining these results into analytical expressions for the HBT radii (see Ap- 
pendix and Section [HI Q ), let us shortly comment on the validity of the transverse saddle 



point approximation introduced in this subsection. We have checked three possible sources 
of deviations of the approximate analytical HBT radii from the results obtained by exact 
numerical integration of the expressions ( |3.6| ): 



(i) The expansion ( ^.7| ) of the Bessel functions K^{a) is asymptotic, but not convergent. In 
our case this is not a problem since we found that in all cases where the approximation 
scheme worked (see below), convergence was reached at low orders p < 3, cf. Section 
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(ii) The factor {ch rjt)^^^^^) in G'„(r) was substituted by exp — ('t- + |) {ch.r]t — 1) , cf. 
Appendix 0. Numerical checks showed that in all the cases considered in Section 
the errors introduced by this approximation are negligible. 

(iii) The crucial approximation in our method turns out to be the substitution of the 
transverse weight function Gn{r) (in the form ([A^)) by the Gaussian ( |3.1CI| ). This 



is good only as long as the quadratic term in the expansion of lnG„(r) around Xn 
dominates the behaviour of the integrand. As we will discuss in Section [VB| , this 
places a strong mathematical constraint on the form of the transversal flow rit{r). 

C. Expressions for the HBT radii 

It is now only a matter of patience to combine Eqs. (|3.3|) , (p.5|) , (pl7|), (|3.13|) , and ( p.l4| ) 
to obtain analytical expressions for all the HBT radii ( p.6|) in terms of the two effective 



Gaussian size parameters ( |3.11| ). This is done in Appendix Since the comparison in 
Section |V] with numerical results will show that we typically must keep several terms from 
the expansion of the Bessel functions, these analytical results are rather complicated and 
don't lend themselves easily to an intuitive interpretation. In particular, different terms of 
the expansion generate different types of dependences, and there is no simple scaling 
law for the HBT radii. 

It is nevertheless instructive to compare our analytical results with simpler ones previ- 
ously proposed by several other authors ||3|J^,^,[7|, |l0tp!3| . Those results were attractive because 
they implied simple scaling laws for the HBT radius parameters as a function of the average 
transverse momentum K± of the boson pair (resp. of its transverse mass m^), and for this 



reason they are still popular with experimentalists (see, e.g., pT|-p^). We will discuss here 
which limits have to be taken in order to recover these results from our expressions, and 
what these limits imply. 

Let us begin with the longitudinal HBT radius Ri. Historically the first expression was 
derived by Makhlin and Sinyukov in a model without transverse flow, via saddle point 
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approximation around r] = r = 0. They found 



T 

= Tq [Makhlin and Sinyukov [10]] , (3.15) 



and this has been used by experimentahsts to fit their data pT|-p3|. For boson emitting 
sources with a non-relativistic hnear transverse flow v{r/R), v <^ c, a. next order correction 
to (|3.15|) was calculated by Chapman et al. who found in the boost-invariant limitj^ ( A77 00 



in their notation) 



i?f = Tq 1 1 + I - + ^ j j [Chapman, Scotto and Heinz [4]] . (3.16) 

In that paper the transverse shift of the saddle point away from f = was neglected, 
which is justified for small transverse flow velocities. An analytical expression without any 
approximation was recently derived for the case of vanishing transverse flow: 

Rf = Tn , \Herrmann and Bertsch [13]] . (3.17) 



In the limit f — > 0, ( ^.l(j| ) reduces to the first two terms from an expansion of ( ^.17] ) for 



m± 3> T. Our expression ( |3.(j| ) for the longitudinal radius in the case of arbitrary transverse 
flow reads (cf. Appendix ^) 

= ^0 /-TT 7^ \ = ^0 5^° ^" f" [our general expression] . (3.18) 



^ For the sake of completeness, we mention another generalization of the Makhlin-Sinyukov for- 
mula, 



T 1 

Rf = Tq ^ [Csorgo and L0rstad 



which was found by saddle point approximation for a source with a finite longitudinal extension 
Ar] which breaks the boost-invariance. In the limit Arj — > 00, this expression reduces to ( 3.15| ). 



The results of Ref. [Q generalize it by non-leading corrections to the saddle point approximation. 
Our investigation here is restricted to boost-invariant models. 
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From this result Eq. (|3.17|) follows immediately in the limit of vanishing transverse flow, 



limr,^=o a{r) = and Eq. ( p.l5|) is the leading term of an expansion of (|3.17|) for large 



T 



■. It was pointed out before that keeping only this lowest order term is an insufficient 



approximation for realistic values of T and ■m±, and that therefore ( |3.15| ) should not be used 
to extract the freeze-out time tq from data. 

In the presence of transverse flow the m^-dependence of Ri is much more complicated 
than either (|3.15| ) or (|3.17|) . The approximate result (p.l6|) for a weak linear transverse 
flow cannot, in this form, be derived from (|3.18| ); one can show that for <^ <^ 1 the 
leading term ~ f ^ agrees in both expressions, but differences occur in higher orders. These 
arise from the transverse shift of the saddle point r which was neglected in Ref. [Q, but 
must be taken into account for stronger transverse flows. This leads to a significant further 
modification of the m^-dependence which will be studied in detail in the following Section. 

Let us now turn to the "side"- and "out" -radii. Here we can again compare^ with the 
expressions obtained by Chapman et al. for a weak linear transverse fiow, expanding 
around r = 0. (Earlier results by Csorgo and L0rstad |^ can be obtained from those of 
Ref. by keeping only the leading terms of this expansion.) In the limit of a boost-invariant 
source — 00 in their notation) and sudden freeze-out at proper time tq (At 0) they 
find in the LCMS frame {(3l = 0) 

J \ [Chapman, Scotto and Heinz [4]] (3.19) 



1 + 



Ri 



R 



Let us compare this with the lowest order [n = 0) contributions from our approximation 
scheme (see Appendix 0): 



^ The results of Akkelin and Sinyukov for systems with strong transverse flow, obtained by 
saddle point approximation, are difficult to compare with since the authors have for mathematical 
convenience modified the geometrical part of the emission function ( |2.3| ) . 
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^si^o) I |-^^^ lowest order contribution] (3.20) 

For a weak linear transverse flow rj = rif{r/R), rjf <^ 1, one may approximate the transverse 
saddle point by xo=0. Then with a little algebra Eqs. (|A^/ |A8| ) can be shown to yield 



n2 



— — — Aq^. Inserting this into ( |3.2U| ) one recovers the first term on the right hand 



A2(0) - A2(0) 

sides of Eqs. ( |3.19| ) except that in the denominator ^ is replaced by Aq = ^ + i. The 
reason for this slight discrepancy is that in Ref. (as in all other previous discussions of 
the emission function ( p.3|) ) the contribution of the factor {chr]t)~^^^ in (3^) to the lengths 
of homogeneity was overlooked. 

These first terms in Eqs. ( p.l9|) arise from the "geometric" parts, (x^) — (x)^ and (y^), 
respectively, in the model- independent expressions (|1.4| ), which for weak linear transverse 
flow coincide to lowest order of our approximation scheme. The second term in the expression 
(|3.19|) for Ro is due to the finite duration of particle emission by the source ( |2.3| ) and 



comes from the term ((t^) — (t)^) in (L4). It is identical with the lowest non- vanishing 



contribution from Eq. (|B6|) . Thus the authors of have correctly identified the lowest order 
contributions to (x^) - (x)^, {y^), and Pl ((t^) - (t)2). (The term -2(3^ {{xt) - (x)(t)) in R^ 
is much smaller and according to ( P3| ) can only be seen if the effect of the transverse flow on 
the saddle point Xj is taken into account.) In the following Section we will show, however, 
that for realistic values of the transverse flow velocity these lowest order contributions are 
unreliable, and a correct description of the i^Tj^-dependence of the HBT radii requires the 
inclusion of higher order corrections from our general expressions in Appendix |B[ 

We already noted that for transverse flows which are not weak or have a non-linear 
dependence on r, the two transverse homogeneity lengths Ao and are no longer equal 
(see also |0). In that case it is no longer guaranteed that the finite duration of particle 
emission yields the dominant contribution to the difference R^ — R^ as in models with weak 
or vanishing transverse flow 0J^j2^ . Since transverse flow is a generic feature of the sources 
generated in heavy-ion collisions [^], its effect must to be taken into account when trying 
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to interpret measured differences between Ro and Rg in terms of the emission time. 



IV. NUMERICAL EVALUATION OF THE CORRELATOR 

The r^-integration for the full 2-particle correlation function C(K, q) as given by Eq. ([1] 
can be done analytically by the same methods as used in Section |111 A| . The resulting 2- 
dimensional transverse Fourier integrals will then be evaluated numerically. The half widths 
of the numerically computed exact correlator can then be compared with the analytical 
expressions for the HBT radii from the previous Section, allowing us to check our approxi- 
mations. 



The 77-integral for the denominator of Eq. (|1.2| ) was already evaluated in Section [III A 



d'xS{x,K) = / rfxrfi/Ki(a(r))e^^(^/^)^^'^^/^e-^'/(2«'). (4.1) 



(27r)3 

For the numerator let us first study the limit qi = 0. Writing t=Tsh.ri, z=Tch.ri and using 
q^=P_Lqo (which is true in the LCMS frame where /3l = 0), we find 

= J dxrfye-*(^°^+^=^)iri(a(r)-zro/3^g„)e^^(^/'^)^'^^'W/^e-^'/(2ij2). (4.2) 

The grdependence of C(K, q) can be obtained by using 

/oo 2 / \ 

d7]chT] e-"^^''-^''^^^ = Ki ( Va2 + 62] (4.3) 

-00 Va2 + 62 V / 

for real values of a and b. This leads to 



j d'^x S{x, K) e 



2rom_L [ a{r) 



dxdy y K, Ua{rf + qfTn e^x^^/^^-^H^)/^ /(.k ) . (4.4) 

With these (exact) expressions the correlation functions C(K, Q's, Q'«=0) and 
C(K, go=Q's=0, g;) reduce to a ratio of 2-dimensional integrals. These have been evalu- 
ated numerically. From the numerical result for the correlator we determine new HBT radii 
i?^(K) by requiring that the function 
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C(K,g,) = l + e-WV, ^ = oJ,s, (4.5) 

reproduces exactly the half point, C(K, q) = 1 + of the numerically computed correlator 
(the other two components qi^j are set equal to zero). In the following Section these values 
Ri(K) will be compared with HBT radii obtained from the expressions (|3.6| ) by either 
numerical or analytical evaluation of the transverse integrals (...)*. 



V. RESULTS 



In this section we give a detailed and quantitative analysis of the HBT radii for the source 
function ( p.l|/ 2.2), both for a linear and a quadratic transverse flow rapidity profile rit{r). 



We will compare the analytic methods developed in Section |T| with exact numerical results. 
We will show that for the linear transverse flow the analytical approximation scheme works 
very well, but only if higher order contributions (non-leading terms in the Bessel function 
expansion and an improved value for the transverse saddle point) are properly taken into 
account. For the quadratic flow the saddle point methods from Section |TTT|are found to fail. 
For a weak quadratic transverse flow we also find an unexpected effect: in the x— y-plane 
the emission region increases with rising transverse momentum of the pair. Although this 
case may be somewhat pathological, it provides a specific counter example to the folklore 
[^,^,^ that transverse collective fiow leads to a reduction of the transverse HBT radii at 
finite transverse pair momentum. Another example with a different velocity profile, but 
similar behaviour for Ro was found in Ref. fTl]. 



A. Linear transverse flow rapidity profile 

All the sources studied in this Section have a transverse geometric (Gaussian) radius 
i? = 3 fm and freeze out along a hyperbola of constant longitudinal proper time tq = 3 fm/c 
at temperature T = 150 MeV. The longitudinal fiow is always boost-invariant, and we study 
the sensitivity to the transverse fiow. 
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Let us begin with a linear transverse flow rapidity profile: 

Vt{r)=Vi^- (5-1) 

For this case the transverse momentum dependence of the HBT radii Rg, Ro, and Ri is shown 
in Fig. |I[ In the figure we compare various approximations: The solid lines are obtained 
by numerical evaluation of the expectation values in the model-independent expressions 
(|1.4| ) in the form ( p.6| ) with our source ( p.l| ). The long-dashed lines show the radii from 
Eq. ( [4 .51 ) which reproduce the half width of the exact correlation function in direction i. The 
short-dashed lines represent our analytical results ( |B2| ) - ([B6|) , evaluated with the effective 
size parameters ( |3.11| ) obtained in Appendix ^ and including all terms up to order p = 3 



from the expansion of the Bessel functions. The dash-dotted lines show the same analytic 



expressions but truncating the expansion (|3.7]) at lowest order p = 0. 



Let us summarize the most important features of Fig. |I|: 



1. For Rs our analytical expression ( P2|) approximates very accurately, even for p = 0, 
the exact value obtained numerically from Eq. (|3.6|) . Thus for Rg the leading order 
expression 

can be used with excellent accuracy. However, it is necessary that the homogeneity 
length \s is evaluated sufficiently closely to the exact transverse saddle point xq as 
described in Appendix |^. Our studies showed that, if used with the lowest order 
estimate (|A^), Eq. (|5.2| ) develops for large transverse flow rapidities r/f > 0.3 a much 
stronger i^^-dependence than the exact "side" -radius. Clearly this renders an ana- 
lytical determination of the transverse flow velocity from the i^T^-dependence of Rg a 
somewhat subtle issue, to say the least. 

The exact numerical value for Rg from Eq. ( ^.61 ) also coincides very accurately with 
Rg from ( [4.5|) . For the source (|2J] ) our model-independent expressions ( |1.4|) thus 
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correctly reproduce not only the curvature of the correlation function at small values 
of qs (which is difficult to access experimentally), but also its half width. This should 
not be too surprising after one checks that in the "side" -direction the source ( |2.1| ) has 
a rather Gaussian shape for all relevant values of K^. 

2. In the "out" -direction the situation is more complicated: there the leading term from 
the Bessel function expansion is seen to yield qualitatively wrong results. After in- 
cluding higher order terms up to p=3, the agreement with the exact numerical results 
improves dramatically, although some deviations remain at large K± for strong trans- 
verse flows. This remaining discrepancy can be traced back to the fact that we use for 
the transverse integration only an approximate saddle point. For the "naive" saddle 
point the disagreement is much worse. This shows that for the "out" -direction and 
strong transverse flow the usefulness of the analytical approach relies crucially on an 
accurate approximation for the transverse saddle point. 

It is important to note that the lowest order approximation (dash-dotted line) com- 
pletely misses the rise of Ro at small values of K±. This is a lifetime effect and arises 
from the time variance given in ([B3|) . The rise was also noted previously in correlation 
functions based on numerical hydrodynamic calculations |26|, but without a detailed 
theoretical analysis of its origin. As seen from Eq. (|B9| ) , the time variance only begins 
to receive contributions at third order of the Bessel function expansion, but the co- 
efficients of these higher order terms are larger than unity. The physics of this effect 



has been explained at the end of Section [IIl7\| . For weak flow (r^f < 0.3) it yields the 



dominant contribution to the difference — Rj. between the "out" and "side" radii. 

For strong transverse flow, that difference remains appreciable even if the contribu- 
tion from the time variance is neglected; it is then due to the difference between the 
homogeneity lengths in the two directions, see Eqs. ([A^/ |A8| ) . This difference can be 
rather accurately estimated by comparing the dash-dotted lines in Fig. |I]a (our lowest 
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order analytical results for Ro which give essentially y (a;^) — (x)^) with Fig. |I|b. For 
?7f = 0.9 one finds, for example, that in the region K± > 500 MeV/c the expressions 
(y^) and (x^) — (x)"^ differ by more than 50%. 

The agreement between the exact HBT radius Ro from the model-independent ex- 
pression ( [1.4|) and the half widths described by R^ is excellent for all values of rjf and 
K±. The reason is again the rather Gaussian shape of the source function also in the 
"out" -direction. 

3. For the longitudinal radius Ri, the leading term from the Bessel function expansion is 
insufficient again (this is even true for vanishing transverse fiow, see ^^), but excellent 
agreement with the exact value from the model-independent expression (|1.4]) is reached 
for p = 3, for all values of K^. For small transverse momenta, however, one sees a 
disagreement with the half width of the exact correlator as given by the parameter 



Rf{K±). On the other hand we checked numerically that Rf{K±) from Eq. ( |1.4| ) 



correctly reproduces the curvature of C{K±, qi) at qi = 0, as it should. In this case the 



disagreement arises from non- Gaussian features of our source (|2.1| ) in the longitudinal 
[rj) direction: for small K± the source decreases at large values of t] much more steeply 
than a Gaussian. However, the discrepancy is small, and in spite of these non-Gaussian 
features of the source the model- independent radii ( |1.4j) , evaluated with the full source 
( p. ID , reproduce the half width of the correlation function with an accuracy of a few 
percent. 

4. All three HBT radii are strongly affected by the transverse flow, even at K± = 0. The 
influence on Ri is relatively weaker than on Rg and Ro- As transverse flow increases, 
the transverse region of homogeneity decreases, and the HBT correlator receives con- 
tributions from a smaller and smaller fraction of the total source. At K± = the 
"out-" and "side-" radii are equal, Rs{K^=0)=Ro{K±=0)] this must be true because 
for = there is nothing to distinguish between the two transverse directions 0. 
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An extraction of the geometric transverse radius R from the data requires to disen- 
tangle geometrical from dynamical effects. Based on Fig. |1| it appears to us that the 
most promising way to achieve this is as follows: one estimates the amount of trans- 
verse flow from the strength of the i^'^-dependence of the "side" -radius Rg. One then 
extrapolates Rs{K±) to = and uses the estimate for the transverse flow rapidity 
to correct it for flow effects, thus extracting the geometric transverse radius R. 

5. From Fig. |l|a it is clear that the extraction of the duration of particle emission, which 
is related to the rate at which Ro{K±) rises for small K±, requires an accurate mea- 
surement of the K^-dependence of the "out"-radius in the domain K± < 100 MeV/c. 
This is not easy, because for very small K± the detection efficiency tends to decrease 
(in a TPC the track density becomes high, in a single arm spectrometer one of the two 
particles tends to miss the detector). It therefore helps to know that Ro agrees with 
Rs at K± = 0; the latter can be extracted much more easily by extrapolation, due to 
its smooth i^j^- dependence. 

Fig. shows that the leading terms from a saddle point approximation produce quanti- 
tatively and qualitatively unreliable results for the m^-dependence of the HBT radii Ri and 
-Ro- Quantitative attempts to extract values for the freeze-out time tq, the transverse flow 
rif and the transverse geometric size R from HBT data must thus be based on a numerical 
evaluation of the model-independent expressions (|1.4| ) rather than on simple m_L-scaling laws 



extracted from insufficient analytical approximations. 

B. Quadratic transverse flow rapidity profile 

A linear transverse flow rapidity profile is perhaps the simplest, but not obviously the 
most realistic assumption. Dynamical studies based on a numerical solution of the hydro- 
dynamic equations [^] produce transverse flow velocity profiles which for small values of r 
look parabolic and for large r saturate at the light velocity; the parabolic rise is zero initially. 
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but builds up gradually as the hydro dynamic pressure generates transverse flow. On the 
other hand, when combining such hydrodynamic calculations with a dynamic model for the 
freeze-out kinetics, taking into account sequential transverse freeze-out of the matter from 
the dilute edge at early times to the center at later times, the resulting flow proflle along the 



at intermediate r and dropping to zero again at the edge of the flretube (because that part 
of the matter freezes out immediately due to its diluteness, before transverse flow is able 
to build up). Numerical kinetic simulations of pion production based on the RQMD code, 
flnally, produce a transverse flow proflle at freeze-out [^] which rises quadratically at small 
r, reaches a maximum and slightly drops again at large r. 

The shape of the transverse flow proflle is thus not at all clear. Investigations of particle 
freeze-out based on hydrodynamic or kinetic models, which result in source functions which 
are only numerically known, and their influence on the shape of the HBT correlation function 
should thus be supplemented by systematic studies of simple model sources with various 
types of simple flow proflles. In this subsection we therefore supplement the results from 
the previous one by a discussion of quadratic transverse flow rapidity proflles: 



The resulting transverse momentum dependence of the HBT radii is shown in Fig. |^. The 
labelling and meaning of the various curves is the same as in the previous subsection. How- 
ever, the curves corresponding to analytical approximations are missing; the reason is that 



rapidity proflle. That scheme relied on a Gaussian saddle point approximation in the trans- 
verse direction; since for small K± the dominant contribution from the transverse flow arises 
from the factor chrjt in the exponent which contributes only at order r^, a Gaussian saddle 
point approximation misses the dominant flow effects completely. An improved analytical 
approximation scheme on the other hand would be much more complicated. 

Comparing Figs. |I| and ^ we see that qualitatively the i^^-dependence of the HBT 



freeze-out hypersurface turns out to be linear again for small r, reaching a maximum 




(5.3) 



the scheme developed in Section |III| and Appendix ^ fails for a quadratic transverse flow 
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radii is similar for both types of transverse flow profile. In all cases the model-independent 
expressions ( |1.4| ) yield HBT radii which give an excellent representation of the half-width 
of the correlator. For a given flow scale rjf, the flow effects on the HBT radii are stronger 
for the quadratic flow ( |5.3| ) than for the linear flow ( |5.1| ). This indicates that emission from 
regions r > i? in the source plays a significant role for the shape of the correlation function. 

The only unusual feature in Fig. |^ is the rise of Ro with K± for a weak quadratic 
transverse flow with rjf = 0.1. It can be traced to a rise of the variance in the "out"- 
direction, (x^) — (x)^, with K±. This is different from the generic decrease of the effective 
region of homogeneity with increasing K± which is seen in all other cases. It seems to be 
due to an accidental coincidence of parameters: In Fig. ^ we show a contour plot of the 
transverse weight function ( p.4|) ; only for small quadratic transverse flow it gets wider in the 
X direction with increasing K_i_. Although exotic, this example should be kept in mind as a 
reminder that transverse flow not always leads to a decrease of the effective source size. 



VI. SUMMARY AND CONCLUSIONS 

The aim of HBT interferometry with hadrons produced in heavy ion collisions is to 
obtain information about the space-time structure of the emitter. Unfortunately, a unique 
reconstruction of the emission function S{x, K) from the 2-particle correlation function 
C(q, K) is impossible, since the frequency g° and the spatial momenta q in the Fourier 
transformation relating the two functions are not independent, as the final state hadrons are 
on mass-shell. However, the class of possible models for the emission function can be strongly 
restricted by a careful multidimensional HBT analysis of the 2-body correlation function. In 
particular, a careful analysis of the dependence of the HBT radii (which describe the width of 
the correlator as a function of the relative momentum q) on the pair momentum K provides 
crucial information on the correlations between emission point x and particle momentum K 
in the source. Such x—K correlations are generated, for example, by the collective expansion 
of the sources generated in heavy ion collisions. Recently heavy ion experiments pT|-p3| have 
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begun to provide first quantitative information on the K-dependence of the HBT correlation 
functions. The hope to obtain access to the collective behaviour in heavy ion collisions by 
analyzing these results has refuelled the theoretical interest in the K-dependence of the HBT 
radii. 

In this paper we reanalyzed this issue with a combination of numerical and analytical 
methods. For a given source model, the numerical methods provide us with exact predictions 
for the transverse momentum dependence of the HBT radii. They thus establish a reliable 
link between various possible models for the emission function and the experimental data 
for the correlation function. The analytical approach provides, in the context of a specific 
source model studied in this paper, a mathematical understanding of the numerical results 
and a bridge between approximate analytical expressions pubhshed previously and the exact 
numerical results obtained here, and thus permits to test their reliability. 

Our numerical studies showed that the model-independent expressions (|L^) for the HBT 
radii not only give the correct curvature of the correlation function C(q, K) for small values 
of q (as was known before Q,]!^), but also reproduce on most cases the width of the correlator 
quite accurately. As long as the measured correlation function has a roughly Gaussian shape 
in q, these expressions provide an accurate and valuable link between theory and experiment 
and permit to relate the (K-dependent) width parameters of the correlation function to the 
(K-dependent) space-time width of the original source function. In this context it should 
be cautioned, however, that the effects of resonance decays, which are numerically known 
to give the correlation function a non-gaussian shape p6| , p8| , remain to be studied in more 
detail. Furthermore, this result means that only the half-widths of the emission function in 
space-time can be reconstructed from the width of the correlation function. Finer details 
of the space-time structure of the emitter (like sharp edges or holes in the center) do not 
affect the width of the correlator, but can only be estimated by a very accurate study of the 
long-range decay of the correlator (e.g. by looking for side maxima at larger values of q). 

The numerical study also enables us to give a detailed assessment of the accuracy of 
previously published approximate analytical results for the HBT radius parameters. This is 
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of particular interest since some of these results have gained great popularity because of their 
simple scaling behaviour as a function of mj_. Here our results cause disappointment: none 
of the so far suggested simple scaling laws is quantitatively reliable, except for very limiting 
situations which are not likely to be established in experiments. To make quantitative use of 
the measured K-dependence of the correlation function and to obtain reasonably accurate 
estimates for the lifetime, transverse geometric size and collective expansion rate of the 
source, a full-blown numerical evaluation of the expressions ( |OD is required, at least. The 
simple scaling laws should not be used - we showed that they yield badly misleading results. 

For the specific source studied in this paper the duration of particle emission is very 
short, of order At{K^) ~ ^JT^ + Rf{K±)/A - tq 0.9 fm/c at K± = 0). This leads to 
a rise of Ro at small K^. The rise stops at K± ~ (when f3± approaches unity), and 
beyond that point Ro and the difference between Ro — Rs decreases again (because At{K^) 
decreases). The experimental verification of the effect of a finite duration of particle emission 
thus requires a high-statistics study of the correlation function at small K±. One might argue 
that realistic models in which the sharp freeze-out along a proper time hyperbola is replaced 
by continuous freeze-out over a longer time period should lead to a larger effect which is 
more easily seen in the data. However, in such models usually the source shrinks as the 
particles evaporate, and particles emitted at later times come from smaller regions near the 
center of the fireball; the net effect on the difference Ro — Rg seems to remain small and 
concentrated at small K± < 100 MeV/c p8|. In the light of these new expectations the 
failure of the experiments so far |]2l|-p^ to provide positive evidence for a finite period of 
particle emission appears less puzzling. We no longer believe that these data show that in 
heavy-ion collisions pion emission occurs "in a flash" . 
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APPENDIX A: DETERMINATION OF THE SADDLE POINT 



Here we give the technical details of our saddle point approximation of expression 
We begin by writing the prefactor as (chr^t)^^"^^) = exp[— (n + |) ln(ch?7i)] and expanding 
the logarithm around chrjt = 1: 

(chr7i)-("+^) = e-^"^^)^^^'"-^) (l + ^^^f + Oir]^)^ . (Al) 

We checked numerically that keeping only the first term in the brackets is sufficient in the 
following. This allows us to combine the exponent with the ch T^t-term from the Boltzmann 
weight factor and to rewrite Gn{r) in the form 

G„(r) ^ exp(n + ^ + t/„(r) - r^/{2R'')) , (A2) 

where 

dn{r) = -Ar,chT]t{r) + B-shT]t{r) (A3) 

r 

with 

^„=n + i + !^, B = (A4) 

The saddle point rj^ = (x^,?/^) of the modified dynamical term exp((i„(r)) alone is found to 
satisfy 

yi = 0, th,^.th„(xf..O) = -^|i^ = |-. (A5) 



Omitting the term + T in the denominator, this reduces to the expression given in [0 
whose shortcomings were discussed in section [V A| . In this case the saddle point is given by 



that point of the source at which the transverse fiow velocity thrji agrees with the transverse 
velocity l3±={K±/m^) of the pion pair. Numerically a gaussian approximation around this 
point leads, however, to unsatisfactory results. Including the effects of the (chr^f)"*^"^^) 
prefactor both in the saddle point and in the gaussian curvature considerably improves the 
approximation. In addition, the saddle point r„ of the full function G'n(r) is further shifted 
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away from by the geometric factor exp(— r^/(2i?^)) in ([A2|) . Due to the symmetry of the 
function Gn{r) this shift will be only in x-direction, i.e. |/„ will stay zero. 

To estimate the position of the exact saddle point r„=(a;„,0) we proceed as follows: we 
expand exp((i„(r)) to second order around an arbitrary point f„=(x„,0). After combining 
the resulting gaussian with the geometric factor exp(— r^/(2_R^)) we try to adjust x„ itera- 
tively to obtain an analytical approximation to the true saddle point x„ of the full weight 
function. 

The quadratic expansion of the dynamical term reads 



dn{xn + Sx,Sy) = - (^Anchrin - Bshr]n) - (A„sh?7„ - B ch.r]n)r]'^Sx 

8x^ by^ 
~2A2(x„) " 1\l(xn) ' 

where rin=r]t{fn), fi'n={dr]t/dx){xn,0), etc., and 

1 

(v4„ shfin- B ch f/„) -^ritiin) 



(A6) 



\ sh ~ 

+ ( A„ ch r7„ - 5 sh r7„) ( — r/t ) +B ^ , ( A7) 

1 



.dy J xl 

{Ar, Shfln- B ch fin) fj" + {An c\lf]n - B sh fjn) fj'n ■ (AJ 



Inserting for x„ the saddle point of the dynamical factor as determined by ([A5|), the 
linear term in ( [A6| ) as well as the first terms on the right hand sides of Eqs. ( [A7| , [A^ ) vanish 
as they should: 

I B^ 1 , 

(A9) 



1 



Al - B^ ivT ■ (AlO) 



At any other point a;„ this no longer true. Inserting the expansion ( [A 61 ) into ( |A^ ) and 
completing the squares we end up with a single gaussian 

lnG„,(r) ^ Ina - ^||^ - , (All) 

whose saddle point x„ is related to the point x„ around which the dynamical factor was 
expanded by 
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•^n -^n 

(A12) 

with 

£n = ^^^n + i?^(n)(A„shr7„ - 5chr7„)C • (A13) 

The radius parameters in (|A1 1|) are given in terms of the geometric radius R and the 
homogeneity lengths (^^,^^) by 



11 111 , 



R^in) R^ Xlixn) ' Rlin) R^ Xl{xn) 
while the constant term is given by 



\nC^ = n+l- (A„chr/„ - Eshr^.) - ^ + . (A15) 



Identifying the expansion point Xn with the true saddle point Xn requires setting £:„=0 and 
amounts to solving the condition ( p.l2|) . The constant J„ then simplifies accordingly. Instead 



of solving ( |3.12| ) numerically we can, however, solve ( [A5| ) analytically for x'^ and then set 



(see Eq. (|ATOD ) 



R^ 



o y^nJ 



This is the saddle point of which would be obtained by expanding (i„(r) around r^. 

We want to stress that this procedure largely corrects for the fact that the radius parameters 
Rj.g{n) should be calculated from the curvature at the full saddle point x„ instead of the 
saddle point x'^ of the modified dynamical factor alone. As shown in Section 0, it is found 
to yield a very good approximation for the HBT radii and their m^-dependence. 

APPENDIX B: CALCULATION OF HBT-RADII 

In this appendix, we give expressions for the HBT-radii ( ^.5p with coefficients calculated 
from the series expansions ( ^.13| / pTT^ ) explicitly up to order p = 3. Using the notational 



shorthand 
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Vm_L/ 



and the results derived in Section IIII BL we obtain 



yp 



,QCnFnRi{n) 



r F 

Z^n=0 '-^n ^ n 



1 3 105 
' 8' 128' 1024' 



and 



Tn 



r F 



1 15 105 \ 
' 8 ' 128 ' ■ ■ ■ ^ • 



The expansion of the "out" -radius -R^ involves three contributions: 

ELo C„ {Rl{n) + Xl) (ELo Cn i^n X^f 



{x') - {xY 



r F 



(Sn=0 F-n 



Elj=od^j'F,F,{Rl{t) + Rlij) 



[Xi 



X 



-2(3^[{xt)-{x){t) 

/5i {{t') - {t) 



{Yln=0 CnFn) 

Xi) 



(T^ c F ? 

Plrmh=od%F,F, 

4 (ELoCnF„)^ ■ 
Here, up to order p = 3 the coefficients are given by 
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-1 

^ 16 



45 
16 



3 
2 
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_9_ 

16 



(B9) 



V / 

It is interesting to note that only the expressions 



and ( p6l) contain coefficients larger 



than 1 (cf. our discussion in Section [III Q ). 
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FIGURES 

FIG. 1. -ftr_L -dependence of the HBT radii Ro (a), Rg (b), and Ri (c), for the emission function 
( |2.3D with parameters T=150 MeV, to=3 fm/c, i?=3 fm. A hnear transverse flow rapidity profile 
r]t{r)=r}{{r / R) was assumed. Curves for different values of r/f are shown. The solid lines are 
calculated numerically from Eq. (|1.4|); the long-dashed lines parametrize the widths R^ of the 
numerically computed correlator C (K_|_ , q) according to Eq. ( [4.5| ) ; the short-dashed lines represent 
our analytical results to order p = 3, while the dash-dotted lines denote the corresponding lowest 
order results. 

FIG. 2. Same as Fig. but for a quadratic transverse flow profile rit{r)=rii{r'^ / R?) . Only the 



exact HBT radii from a numerical integration of the mo del- independent expressions (1^) (solid 
lines) and the width parameters i?? of the numerically computed correlator C(K_|_,g) according to 
( [4. 51) (dashed lines) are shown. 



FIG. 3. Contour plots for the emission function (|3.4| ) in the transverse x—y plane, with pa- 
rameters T=150 MeV, ro=3 fm/c, R=2> fm, for a quadratic transverse flow profile r]t{r)=ri^{r'^ / R^) . 
From center to edge the lines correspond to 90% . . . 10% of the peak value (in steps of 10%). The 
two left diagrams show that for a weak quadratic transverse flow with r^f = 0.1) the emission region 
actually increases in the x-direction with increasing transverse momentum. The two right diagrams 
(for r]{ = 0.3) show the generic decrease of the source in both x and y directions with increasing 
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